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Abstract. A lattice system is derived which amounts to a higher-rank analogue of the Q3 
equation, the latter being an integrable partial difference equation which has appeared in the 
ABS list of multidimensionally consistent quadrilateral lattice equations. By construction this 
new system incorporates various lattice equations of Boussinesq type that were discovered many 
years ago. A corresponding Lax representation is derived. 



1. Introduction 

In recent years the integrability of partial difference equations (PAEs) on the two-dimensional 
space-time lattice has become a subject of considerable interest. The early examples of such 
equations date back to the 1970s [B and early 1980s [TOJ E2J gl], but it is only in the last 
decade that their study has moved to the centre-ground of interest in integrable systems. The 
notion of multidimensional consistency, cf. [23] [8], as a key integrability characteristic, has 
formed the trigger of a wealth of activity in the subject, leading among others to the celebrated 
classification result by Adler, Bobenko and Suris (ABS), [2], for scalar quadrilateral PAEs. For 
the equations in the ABS list (which also comprise some of the quadrilateral lattice equations 
already found in the previous era) many results are known, in particular on special solutions, 
cf . [U (5J El [I9j [20"1 [7] . A classification result was also formulated recently for three-dimensional 
scalar octahedral equations, i.e. equations of lattice KP type, [3j, but is markedly still absent 
for multicomponent and higher-order lattice equations. In particular it seems that classification 
results for the "next" class of equations beyond the scalar quadrilateral case, i.e. those of 
Boussinesq type, requires new ideas and new technical concepts for its treatment. 

As far as I am aware, a first example of a lattice Boussinesq (BSQ) equation appeared in |10j . 
this equation being a straight dimensional reduction of the Hirota bilinear KP equation, |14j . 
it is not clear to me that it is rich enough to allow for all appropriate (including intermediate) 
continuum limits. It is a special case of a more general lattice BSQ equation which appeared 
in [21] where the general class of higher-rank lattice equations of what we coined "of Gel'fand- 
Dikii type" was considered. (This class, labelled by roots of unity oj = exp(2iri/N), comprises 
the KdV type lattice systems for N = 2 and BSQ type systems for N = 3). The lattice BSQ 
equation, as given in [21], can be presented in the form of the following 9-point equation on the 
two-dimensional lattice: 

p3 _ q3 p3 _ q3 



P — q + n n+ i im+ i — U n+ 2,m P — q + Un,m+2 ~ '"n+l.m+l 

= (p-q + u n+ i im+2 - u n+2 , m +i)(2p + q + u„ im+ i - «„+2,m+2) 

-{p-q + ^n,m+l ^n+l,m Un,m ^n+2,m+l ) • (1-1) 
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in which u = u nj7n denotes the dependent variable of the lattice points labelled by (n, m) € Z 2 . 
In (11. ip the p, q are continuous lattice parameters associated with the grid size in the directions 
of the lattice given by the independent variables n and m respectively. For the sake of clarity I 
prefer to use a notation with lattice shifts denoted by 



U — Un,m 1 ? U — U ra _|_i jm , U — Un,m 

in terms of which we have also 



H- ^ U — Un^m+i 



U = U n+ i tTn+ i . U = n n +2,m+l , U = U n +l, m +2 , U = U n +2,m+2 

Thus, the lattice BSQ equation (jl.ip is associated with the following stencil 



u 



u 



u u u 

In the same paper |21j we gave a lattice version of the modified BSQ equation, whilst in [17] 
the Schwarzian lattice BSQ was presented. All three equations are specialisations of a more 
general lattice BSQ type equation, which was presented in \17\ [T8] and which can be cast into 
the following form: 

P a Pb(uu + uu) — Q a Qb(uu + uu) — (j> 3 — q s )(uu + uu) 



P a Pb{uU + UU) — QaQb(uU + UU) — (j) 3 — q 3 )(uU + Uu) 

_ {QgU - Q b U) (PgU - P b u) (QgP b U - PgQbU) 

(P a u - P b u) (Q a u - Q b u) (Q a P b u - P a Q b u) 
In eq. (I1.2p the P a , P b , Q a and Q b are given by 

Pa 



(1.2) 



Vp 1 



6 3 



Qa 



Qb 



b 3 



where a and b are some additional (continuous) parameters. If b = a (11.21) reduces to the lattice 
Schwarzian BSQ equation, whilst in the limit b — > oo, bu — > v, we obtain the modified BSQ for 
the variable v. In the double limit a, b — > oo we recover (jl.lh after an appropriate scaling. Thus, 
(|1.2j) can be viewed as an equation interpolating between these various cases, in the same way as 
the so-called NQC equation of [22] interpolates between various lattice equations of KdV type. 
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The mentioned lattice BSQ equations can also be written in multi-component form, and 
as was shown in [18] the modified BSQ can be written as 2-component system of PAEs on 
an elementary quadrilateral, and as such multidimensional consistency in the sense of [23[ [8] 
can be easily shown for that system. In [27} [25] eq. (jl.lj) was written as a three-component 
system on the elementary quadrilateral, whilst the multidimensional consistency for the 9-point 
scheme of the scalar form was studied in [27]. Some further multicomponent BSQ type systems, 
generalizing the 3-component version of the lattice BSQ, were given in a recent paper jllj . We 
mention also that explicit soliton solutions of the lattice BSQ equation were constructed in 
[121 fl6] , although it must be said that soliton solutions can also quite easily be inferred from 
the direct linearization scheme presented in [21} [27] . 

In the present paper I construct what I believe is a novel lattice system which incorporates the 
previously known BSQ type equations, including the modified and Schwarzian BSQ equations 
and (|1.2p . The construction of this system is based on the structure of the Q3 solutions, which 
is a pivotal equation in the ABS list, [2]. In fact it was shown in [19], cf. also [61 [20], that 
solutions of Q3 can be written as a linear combination of four solutions of the NQC equation 
of [22], which is the analogue of (jl.2p in the quadrilateral scalar case. Here I perform a similar 
construction for the higher-rank case, leading in this case to a non-autonomous lattice system 
which I believe constitutes the proper analogue of the Q3 equation. Since, at the present moment, 
even a preliminary classification of higher-rank or multicomponent lattice equations is lacking, I 
believe these results could form a first step towards an insight into the structures behind lattice 
equations in the higher-rank case. In section 2 I will present the basic construction of solutions 
for Q3, which motivates the generalization to the BSQ case. The latter case is markedly more 
complicated, and in section 3 the system of basic constitutive relations is derived, from which in 
section 4 I derive the relations forming the system. Multidimensional consistency is built into the 
system by the entire construction, but the integrability is made explicit by presenting in section 
5 a Lax representation. The whole construction can, in principle, be generalized naturally to 
the whole Gel'fand-Dikii hierarchy, of which the BSQ system represents the N = 3 case. 



2.1. Infinite matrix scheme. The starting point of the construction is the following integral 



over a yet unspecified set of contours or arcs Tj in the complex plane of a variable k, and in 
which for N > 2 a positive integer, uj = exp(27ri/iV") is the N th root of unity. The factors p k , 
aw are discrete exponential functions of k, k' respectively, given by 



p k (n,m) = (p + k) n (q + k) m p k (Q,0) , a k ,{n, m) = (p - k'y n {q - k') m M0, 0) , (2.2) 



and the infinite-component vector c k = (k J )j^z denotes the vector of a basis of monomials in 
the variable k. The integration measures d\j(k) remain unspecified, but we assume that basic 
operations (such as differentiations w.r.t. parameters or applying shifts in the variables n and 
to) commute with the integrations. 



2. Constitutive system 
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In accordance with the notation introduced in section 1, we denote shifts over one unit in the 
variables n, and m respectively by the operations ~ and ~, which implies for the factors p k , cr k > 

Pk = Pk{n + l,m) = (p + k)p k , p k = p k (n,m + l) = (q + k)p k , 

and 

0V = cr v {n + l,m) = (p — k'^ay , d k > = a k ,(n,m + I) = (q - k')~ l a k i . 
These relations imply the following linear relations for the matrix C 

C(p- t A) = (p + A)C , C(q- t A) = (q + A)C, (2.3) 

where we have introduced the matrices A and t A which are defined by their actions on the 
vector c, and on its transposed vector, as follows: 

A c k = k c , c\, A = k' c k , . 

The following ingredients determine the structure: 

i) A Cauchy matrix Q obeying the relation CIA + *Afi = O , where O is a rank 1 projection 
matrix, obeying O 2 = O. 

ii) A matrix U obeying the relation U = C — U ftC 
in) A vector u k defined by u k = p k (c k — U Ct c k ) 



In terms of these objects the following sets of relations involving the shift can be derived: 



U 



U(p- *A) 

N-l 

H(cjip- *a) 



(p + A)U -UOU , 

N-l 

H(oJp + A) 

i=i 



N-2 



j=0 lk=l 



n(«* P -*A) 



O 



N-l 

n (u k p+A) 

k=j+2 



N-l 



U(-*A) N 



A N U - U (" t A) j OA N ~ 1 - j U , 

j=0 

(see [21] )■ By virtue of the covariance of the dynamics in terms of the variables n and m, 
similar relations to (j2.4ajl and (I2,4b[) with the shift ~ replaced by ^ while replacing p by q. The 
associated linear problems (Lax pairs) are derived in terms of the object u k , for which we have 
the following set of constitutive relation, |21j . 



(2.4a) 
U , 
(2.4b) 
(2.4c) 



u k 



(k 



N 



p) N )u k 



(p + A)u k -UOu k , 

JV-l 

H(u>p+a) 

i=i 



N-2 

u k -UY, 

3=0 



k=l 



o 



N-l 

n (w*p+a) 

k=j+2 



N-l 



(2.5a) 
u k , 
(2.5b) 
(2.5c) 



k N u k = A N u k -Uj2 (" WOA"- 1 '* u k , 

3=0 

and similar relations to (|3.2ap and (j3.2bj) with the shift ~ replaced by ^ while replacing p by q. 
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For any fixed value of N, these abstract equations form an infinite set of recurrence relations 
denning the dynamics in terms of the independent variables n, m on the objects U and 
taking values in an abstract vector space V, and an adjoint vector *e in its dual V*. To give a 
concrete realisation, we can choose a fixed vector e and use the matrices A and *A to define a 
natural gradation in this vector space, In terms these, we realize quantities such as C and U as 
infinite-dimensional matrices. In fact, setting 

Ui 4 := t e t K i UK i e , i,jeZ, 

we obtain infinite, i.e., ZxZ matrices, on which A and *A act as index-raising operators acting 
from the left and from the right. In this realization the matrix O = e e is the projector defined 
by these vectors, and this defines the "central" entry in the space of Z x Z matrices. We will 
define concrete objects in terms of the matrix U, where we allude to this realization. 

2.2. Basic objects and their relations: KdV case (N = 2). Thus, for the case N = 2 the 
relations (12.41) take the form: 



whereas (12.51) become 



!7(p-*A) = (p + A)U-UOU, (2.6a) 
U(p+ t A) = (p-A)U + UOU, (2.6b) 
U *A 2 = A 2 U - U (O A - *A O) U , (2.6c) 

Uk = (p + A)u k - UO u k , (2.7a) 
(p 2 - k 2 )u k = (p- A)u k + UOu k , (2.7b) 
k 2 u k = A 2 u k — U (OA — *A O) u k , (2.7c) 

Since (I2.6bj) can be obtained from (|2.6aj) by matrux transposition, the system admits solutions 
that are symmetric with regard to this operation, i.e., U = U . 

We shall briefly show how the Q3 equation, [2], arises in a simple manner, presenting the 
essence of the derivation in [15] . In fact, we need to introduce the following quantities: 

v a :=l- t e(a + A)" 1 Ue = l- t eU(a+ *A) -1 e , (2.8a) 

s a :=a- *e (a + A) -1 U t A t e = a — t eAU(a+ *A)^ e . (2.8b) 

as well as the quantity no := C/o,o and 

s afi := *e (a + A)" 1 U (b + *A)- 1 e = s fe>a (2.9) 

Here a and 6 are arbitrary scalar parameters, and we have suppressed the multiplication by a 
unit matrix in expressions of the type a + A := al + A. The following relations follow from 
(pHaD and (l2~6b|) : 

s a = (p + iio)«a - (P - a)t> a , (2.10a) 
s« = (p + a)?; a - (p-«o)«a (2.10b) 

as well as the relation 

1 + (p - a)s a ,b ~ (P + &)Sa,f> = Vaffe (2.11) 

and using the symmetry of s a) 6 we have another relation of this type with a and b interchanged. 
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Introducing the plane wave factors 



p a (n,m) 



p + a 
p — a 



q + a 



/9 a (0,0) , Pb(n,m) 



p + b\ n ( q + b 



p — b J \q — b 



Pa (0,0) 



we can now introduce the two-component vectors: 



1/2 
Pa V a 
1/2 
P-a v ~a 



1/2 
Pa S a 
1/2 
P-a s ~a 



and the quantities: 



u :- 



P l j\p^ a 



CS—n.h 




where S a h =: s 



U := vl Av b . 

l/(a 4- b), and where A is the matrix: 

A B 
C D 

The relations ()2.10|) lead to the vector relations 

Sa = (P + U Q )v a - P a V a 

{p - u )v a 



(2.12) 

(2.13a) 
(2.13b) 



P a v a 



(2.14a) 
(2.14b) 



and similar relations for the dynamics in the variable m, with the usual interchangements and 
changing P a into Q a , where 



Pa = V? 



Qa = \/q 



and we will also introduce the quantities P = P a Pb , Q = QaQb ■ 

The relation (|2.1ip and its counterpart interchanging a and b, now become 

P a u - P b u = v l a Av h P h u - P a u = v l a Av b , 

where it should be noted that the matrix A need not be symmetric. 

The construction of the Q3 equations from these data now follows the following steps: 
i) Biquadratic identity: Consider the following 2x2 matrix 



(2.15) 



E, 



a,b 



v l a Av b , v^Asb 
slAvh , stAs h 



(2.16) 



The determinant of this matrix can be computed in two different ways; on the one had we have, 
using (|2.14p the sequel: 



det(£ aj6 



Ph 



s\Av b 
U 

P a u - P b u 



v a Av b 

SnAVh 



PaPb 



v a Av b 
viAvh 



v a Av b 
v a Av b 



P b u - P a u 
U 



P 



UU - P{u 2 + u 2 ) + (2p 2 - a 2 - b 2 )uu 



A HIGHER-RANK VERSION OF THE Q3 EQUATION 



On the other hand, using a special version of the Weinstein-Oranszajn formula, the determinant 
equals 

det(E ajb ) = det(A) det(v a , s a ) det(v b , s b ) = a/3 det(A) , 
where a := det(t> a , s a ) and f3 := det(v b , s b ) are constants, as can be seen from the chain: 
a = det(v a , s a ) = -P a det(v a , v a ) = det(v a , P a v a ) = det(v a , s a ) 



a 



using once again the vector relations (|2.14p . and a similarly for f3, and the same argument holds 
for the shift in the variable m. 

Thus, equating both sides of the evaluation of det(E ab ) we obtain the factorization of the 
biquadratic 



UU = P(u 2 + u 2 ) - (2p 2 



b 2 )uu + ^p- det(A) , 



(2.17) 



ii) Q3 equation: To obtain the main equation of interest, Q3, for the variable u we proceed as 
follows. Consider the relations (|2.15|) and insert them into the following determinant 



v l a Av b , v l a Av b 
v a Av b , v a Av b 



P b u - P a u 

QaU ~ Q b U 



QbU ~ Q a U 
P a U - P b U 



= P(uu + uu) — Q(uu + uu) — (p 2 — q 2 )(uu + uu) 
On the other hand we have that the starting determinant equals: 



1 

Qb 



v a Av b , v a A [P b v b - (p-q + u - u )v b ] 
v a Av b , v a A [P b v b - (p-q + u - u )v b ] 



p-q + u -u 
PbQb 

p-q + UQ-up 

PbQbQa 



V a As a , V a Av b 



v a As b 



v a Av b 
v a As a 



-P a v a + (p + q + u - u )v a ) As 



v t a Av b 

-P a v a + (p + q + u - u )v 



(p-q + UQ- u )(p + q + u -u ) 



v a As b , v a Av b 
s\As b , s a Av b 



PaPbQaQb 

Here we have made use of the intermediate relations 

P a V a ~ QaVa = (p-q + U0- U )V 
P a V a + QaVa = {p + q + U - Uq)v 

as well as of the following equation for uq 



p 2 — q 2 
PQ 



a(3det(A) . 



(p - q + Uq - U Q )(p + q + U - U ) 



p 



(2.18) 



which is the HI equation (lattice potential KdV equation) in the terminology of the paper [2]. 
Thus, we obtain in the case N = 2 the Q3 equation from a simple structure of determinants. 
We will now follow the same path for the case N = 3, leading to a higher-rank case, but we will 
see that the structure is considerably more involved in this case. 
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3. Higher rank case 

3.1. Basic objects and their relations: BSQ case (N = 3). In the case N = 3 the set of 

relations 

Thus, for the case N = 3 the relations (12.41) take the form: 



I7(p-*A) = (p + A)U-UOU , (3.1a) 
U{cop- l A){uj 2 p- *A) = (wp + A)(w 2 p + A) 17 + J7 [pO -(OA- *A O)] U , (3.1b) 
t/(-'A) 3 = A 3 £/-£/ (OA 2 - *AOA+ *A 2 0) ?7 , (3.1c) 

and similar relations to (|2.4a|) and (|2.4bp with the shift ~ replaced by ~ while replacing p by q, 
as well as 

Uk = (p + A)u k -U Ou k , (3.2a) 
(p 3 + /c 3 )u fc = (wp + A){uo 2 p + A)S fc + U [pO - (O A - *A O)] S fc , (3.2b) 
k 3 u k = A? u k -U (OA 2 - *A O A + *A 2 O) w fe , (3.2c) 

(and similar relations to (|3.2ap and (|3.2bp with the shift ~ replaced by while replacing p by 
q). For the sake of obtaining from the constitutive relations (|3.ip closed- form equations, we 
introduce the following objects: 

Va ■- l _ t e (a + A)~ l Ue , w b := 1 + 'e U (6 - 'A)" 1 e , (3.3a) 
s a := a - *e (a + A) -1 U l Ae , i 6 := -b + *e A t/ (6 - 'A)^ 1 e , (3.3b) 
r a ■= a 2 - *e (a + A)" 1 ?7 *A 2 e , z b := fr 2 + t eA 2 U(b- 'A)" 1 e , (3.3c) 

but in particular the following object: 

s a , 6 := *e (a + A)" 1 [/ (-6 + 'A^e (3.4) 

For these objects we have the following relations 

1 + (p - a)s atb - (p - b)s a>b = v a w b , (3.5a) 

(p + a + b) + (p 2 +pb + b 2 )s a M - (p 2 + ap + a 2 )s aib = s a u7 b - u i 6 + pv a w b , (3.5b) 



as well as 



uo,i + «i,o = P^o - P^o + «o«o , (3.6a) 
- p 2 +pq + q 2 - (p + g)(«o - u ) - uqvlq + ui o + «o.l , (3.6b) 



p - g + u - n 



and 



and 



s a = (p + u )v a - (p- a)v a , (3.7a) 

t b = (p - b)w b - (p - u )w b (3.7b) 

r a = ps a - (p - a)s a + u a M ,i , (3.8a) 

z b = (p - 6)f b - pt b + ui fi w b , (3.8b) 
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and 

r a = (P 2 + ap + a 2 )v a -(p-u )s a -(p(p-u ) + u lfi )v a , (3.9a) 

Zb = (p 2 + pb + b 2 )w b + (p + u )t b - (p(p + u ) + u 0: i) w b . (3.9b) 

All relations (|3.5p - (|3.9p also hold for their counterparts obtained by replacing the shift ~ by 
the shift *~ whilst replacing the parameter p by q. We mention that from the relations (|3.6p by 
eliminating the quantities uo,i and ui o we obtain immediately the lattice BSQ equation (jl.ip 
in the form: 

po _ go p° — q° 



p- q + u -u p- q + u -u 

= (P ~ q + u -u )(2p + q + u -UQ) - (p - q + u - u )(2p + q + u -u ) . (3.10) 
3.2. Basic vectorial relations. Introducing the 3-component vectors 

PaVa \ ( CF b W h 

PcuaVcua , W b := [ <Tub™0jb (3.11) 

. Puj 2 aVta 2 a J \ C^b^^b 

where 



K p-aJ \q-aj \ P h J \ Q b 

and similarly the vectors s a , t b , r a and z b , we obtain the sets of vectorial equations: 

s a = (p + u )v a - P a v a (3.13a) 

t h = P b w b - (p - u )w b (3.13b) 



as well as 



and 



r a = ps a ~ P a s a + u ,iv a (3.14a) 
z b = P b t b -pt b + u lfi w b . (3.14b) 



r a = P a v a - (p - u Q )s a - {p{p - uq) + uifi) v a (3.15a) 
z b = P b w b + (p + u )t b - {p(p + uq) + u 0) i)w b . (3.15b) 
The main object of interest is the following quantity 

2 

u := 

i,j=0 

for which, as a consequence of (|3.5ap and (|3.5bj) respectively, we have the following basic relations 



2 1 
^2 A i jPui a Oui b S(uJ l 'a,ujib) where S(a,b) = s a:b — , (3.16) 



P a u-P b u = v^Awf, , (3.17a) 

P b u-P a u = slAwb-vlAtb+pv^Awi, . (3.17b) 

It is in terms of the quantity u that we will obtain the higher-rank (BSQ) analogue of the Q3 
equation. We will now show how the equation for this object is constructed. 
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Introducing the quantities: 



u 


= v t a Aw b , 


V 


p = 


= PaPb j Q = 


QaQ 


Pu 


-Qu- (p 3 - 


q 3 )u 


Pu 


-Qu- (p 3 - 


q 3 )u 


Pu 


-Qu- (p 3 - 


q 3 )u 



At b , W := s l a Aw b , Z := s l a At b 
have the following relations 

-(p- q + uo-uo) (p + q-u )U + W , 
(p-q + u -uo) (p + q + u )U-V , 



Z-(p + u )V + (q-u )W 



+ 



p{jp + uq) + q(q - u Q ) + «i o + Uq i 



U 



(3.18) 

(3.19a) 
(3.19b) 

(3.19c) 



which we will think of as Miura type relations. Similar relations appeared in |19j connecting 
solutions of Q3 to solutions of HI (the potential KdV equation), which is the N = 2 analogue 
of the lattice BSQ equation ()3.10p . 

Several additional relations can be derived by combining the relations ()3.13|) , (|3.14p and ()3. 15|) , 
namely from (|3. 13[) by elimination of s respectively t we get: 



PaVa ~ QaV a 

P b w b - Q b w b 



ip-q + uo- u )v a , 
ip-q + uo- u )w b . 



Similarly, by combining (|3.14p in two different lattice directions we find 

PaVa-QaVa = (P ~ Q + ^0 ~ «o) S a + ip + q - Uq)v a 

P b w b -Q b w b = -(p-q + uo-uo) t b - ip + q + u )w b 



Furthermore, by combining (|3.15j) in two different lattice directions we find 

PaSa-QaSa = {p ~ q)s a + [(p - q + U - U )u ~ pu + qu ]v a 



Pht 



Qbtb 



ip - q)t b + ip - q + u - uq)u - pu + qu 



w b 



(3.20a) 
(3.20b) 



(3.21a) 
(3.21b) 

(3.22a) 
(3.22b) 



We will now proceed by using these constitutive relations to derive the lattice equations, involv- 
ing primarily three objects, u, U and uq. 



4. Lattice BSQ analogue of Q3 

The main object in the construction, following the relations obtained in the previous section, 
is the following 3x3 matrix 



E a ,b '■- 



v\Aw b 
siAw b 
r{ A w b 



v* a At b 
slAt b 
r* At b 




(4.1) 
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The determinant of this matrix can be expressed in two different ways. On the one hand, by 
applying a special case of the well-known Weinstein-Aronszajn formula, we have 



det(E, 



a.b) 



det(A) det(v a , s a ,r a ) det(w b , t b , z b ) 

Pa \ ( Qa 



Ph 



Qi 



a/3 det(A) 



(4.2) 



in which a and j3 are constants (w.r.t. the variables n and m) but will depend on the parameters 
a and b respectively. In fact, by applying (|3.14ap and subsequently (|3.13ap and (|3.15ap we have 

det(v a , s a ,r a ) = -P a det(v a ,s a , s a ) = P a 2 det(v a , v a , s a ) = P a det(r a , v a , s a ) 

det(« a , s a , r a ) = P a det(t> a , s a , r a ) , 

similarly we find in the other lattice direction det(u a , s a , r a ) = Q a det(w a , s a , r a ) . Thus, we 
have 

det(v a , s a ,r a ) = P™Q™ det(««, s a , r o )| n=m=0 =: PaQ^a , 

setting the value of this determinant at n = m = equal to a. By a similar argument we find 
that determinant det(w b , t&, z&) can be found to be equal to 

det(w b ,t b ,z b ) = P b - n Q- m det(w bl t b ,z b )\ n=m=0 =: P b n Q b m p , 

setting the value of the latter determinant at n = m = equal to f3. 

On the other hand, using the relations f|3. 13j) (|3.14p we can rewrite the determinant of E ab 

as 

U , v^Awb , v^Awb 
v^Awb , v^Awh , U 



det(E aib ) = -PQ 



(4.3) 



vlAw b , U , v l a Aw b 

(the under-accents ^ and ^ denoting the shifts in the opposite directions to the shifts ~ and 
where the entries can be evaluated using the relations 

Q a P b u - P a Q b u 



v a Aw b 
Pa-vi A w b 
PbV a A w b 
PbQbvl A w b 
P a Q b vlAw b 



p - q + n - iio 

(2p + q + u - u )U + 

(2p + q + u ~ u )U - 



Pu — Qu — (p 3 — q 3 )u 
p - q + UQ -uq 

Pu — Qu — (p 3 — q 3 )u 
p-q + uo-uo 



PU -QU 



P 



+ UQ-UQ 

3 „3 



+ Uq - Uq 



+ Pu + Qu- (p 3 + q 3 - a 3 - b 3 )u 



U - Pu + Qu + (p 3 - q 3 )u , 



(4.4a) 
(4.4b) 
(4.4c) 
(4.4d) 
(4.4e) 



as well as, where relevant, their counterparts obtained by interchanging the shifts and 
and simultaneously interchanging p and q. Eqs. (|4.4p can be deduced from (|3.19p and (|3.13|) 
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through the definitions (|3. 18j) . using also (I3.17P and the relations (|3.20p - (l3.22p . In addition to 
these relations, we also obtain the following important relation 

PbQaU ^QbPaU = p a Q aU + P b Q b u - Q a P b U - P a QbU , (4.5) 

p - q + u - u 

accounting for the internal consistency of the relations (|4.4|) . Eq. (|4.5|) can be used to eliminate 
the canonical combination 7 : = p — q + UQ — uq from all relations (14.4p in favour of the variables 
U and u. Let us now formulate the main statement: 

Proposition: Eq. (|4.2p , where det{E a ^) can be rewritten in various forms such as the one 
exhibited in (|4.3p constitutes the N = 3 higher rank generalization of the Q3 equation, where the 
entries of the determinant can be written in terms of the "conjugate " variables u and U using the 
relations in f|4.4j) (and their counterparts, where relevant, interchanging p and q and~ and~ ). 
The system is supplemented by a set of linear ordinary difference equations, in each of the lattice 
directions, for U obtained by eliminating between (|4.4bp and (|4.4c|) (and their counterparts in 
the other lattice direction) the left-hand sides. In all expressions the canonical quantity 7 is 
obtained from (|4.5p in terms of u and U. 

. We will refrain from writing out the resulting system in explicit form, which may not be 
insightful. The presentation of the system may depend, in fact, on the choice of the combination 
of dependent variables considered as the prime variables. For instance, U can actually be solved 
from a cubic equation which is obtained from det(i£ &) by using the linear difference equations 
for U (alluded to in the proposition), to eliminate U and U from (|4.3[) . We intend, in a future 
publication, to address these options in more detail in order to find a convenient representation 
of this system which is incorporated in (|4.2j) and the relations (|4.4j) , (14. 5p . 



5. Lax representation 

The Lax representations can be obtained in a systematic way from the relations (|3.2p . i.e. 
in terms of quantities defined on the basis of the vectors u k yielding the objects obeying the 
relevant linear systems of equations. The first step is to define the quantities 

4 J) := l eh?u k , u k (a); = l e(a + A) -1 u k , (5.1) 
for which we can derive directly from (13. 2p the following relations 

4 0) = (p - 5 o)4 0) + 4 1} » ( 5 - 2a ) 

Sfe(a) = {p - a)u k {a) + v a uf\ (5.2b) 

(p 3 + k 3 )u k (a) = {p 2 + ap + a 2 )u k {a)-{pv a + s a )u[ 0) +v a u k 1] , (5.2c) 

and a set of similar relations with p replaced by q and ~ by ~. The next step is to introduce the 
3-component vector 

/ p a u k {a) \ 

u k (a) = p^u^a) , (5.3) 
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with the plane- wave factors p a as given in (|3.12p . With the definitions of the 3- vectors (13.111) 
we now obtain from (|5 . 2f) the following relations 



u k (a) 



P a U k (a) + VaU^ 



(5.4a) 
(5.4b) 



(p 3 + k 3 )u k (a) = P a u k (a) - (pv a + s a )ul> + v a u k 

with similar ones involving the other shift of the lattice. 
The third step is to introduce the scalar object 

u k (a,b) := w\ A 1 u k (a) , 

where A t is the transpose of the 3x3 matrix A used in (|3. 18[) . Multiplying the relations ([5 
from the left by w l b and w b respectively, we obtain the set 

w\ A 1 u k (a) = P a u k (a, b) + (w\ A 1 v a ) u[ 0) , 



u k (a,b) 
(p 3 + k 3 )u k {a,b) 
(p 3 + k 3 )w t b A t u k (a) 



P> 6 A V>) + t^f , 

P a w\ A 1 u k {a) - (pU + W)uf + Uu[ l) , 

P a u k (a, b) - {pwl A t v a + w\ A* s a ) u[ Q) + {w\ A 1 v a )u k 1] 



Upon elimination from this set the quantities w t b A t u k (a) and w t b A t u k (a) we obtain the fol- 
lowing relations 



(a 3 + k 3 )u k (a,b) 
(a 3 + k 3 )u k (a, b) 



PaiwlA'v^uf 



(pU + W)u[ Q) + Uu 



r(i) 



(5.5a) 



(p d + k d )Uu k 
Introducing, finally, the 3-component vector 

(f> k (a,b) : 



P a {pwl A 1 v a + w\ A 1 s a ) uf + P a {w\ A 1 v a )u [l) 



k ■ 

(5.5b) 



4 0) 



\ u k (a,b) 

we now obtain the following set of linear equations (Lax pair) for cf) k (a,b): 

4> k (a,b) = L k (a,b)4> k (a,b) , 4> k {a,b) = M k (a,b) 4> k (a,b) 

with the Lax matrices L k and M k . L k is given by 

(p - uo)U , U 

(p-u ){pU + W)-P a E , pU + W 



L k (a,b) 



U 





+ k 3 



, -^P a (TW-UA) , p a r 



(5.6a) 



(5.6b) 



in which the espression * stands for 



UU + 



1 



(p 3 - a 3 ){UU - ET) + P a (p- u )(TW - UA) 



a 3 + A: 3 I 

i.e., such that the determinant is given by det(L k (a,b)) = (p 3 + k 3 )U/U . Here we have used 
the following abbreviations: 



E := wlA t v a 



T'-wlA^a , 



A := wlA^a , 



li 
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which can be identified with the objects given in (|4.4p leading to 

E = P a u-P b u , P a T = (p + u )U -W . P a A = pW + u 0A U - R , (5.7) 

where we have set R := r^Aw^ . The latter quantity can be expressed in terms of more familiar 
objects, such as 

R=(p 3 - a 3 )u -Pu-(p- u )W - {p{p - n ) + Si, ) U = pW + u ,iU , 

The matrix is given by similar expressions making the usual replacements, and the 
compatibility condition of the Lax pair, i.e., the Lax equations L^M^ = M^L^, can be explicitly 
verified entry-by-entry and leads to a set of relations compatible to the relations (|4.4p . Once 
again, for the system the question arises as to what is the best realisation of the system, in order 
to obtain an effective Lax pair which could be used, e.g. for the construction of wide classes of 
solutions of the system. 

6. Conclusions 

In this paper we have derived a new integrable lattice system of higher rank which incorporates 
the special subcases of known BSQ type solutions as were presented in [51] . In fact, by choosing 
the arbitrary matrix A to have a special form, such that the inner products of the type x t Ay 
factorize, it can be shown that the system reduces again to equations of the form (|1.2|) . A new 
feature is that whenever det(A) ^ the analogous term associated with the 5 2 -parameter (in 
the form of Q3 as given in the ABS list of [2]) becomes non- autonomous. This is a feature that 
we expect to persist in the higher members of the Gel'fand-Dikii hierarchy as well, and only if 
N = 2 this "deformation term" of the equation is independent of the variables n and m. We 
stopped short of writing the lattice system in fully explicit form for reasons explained earlier, 
and finding a convenient explicit realisation would be a next step. Many questions are obviously 
worth considering, such as the construction of elliptic solutions along the line of [20] (which 
I expect would lead to basically the same system of equations), the question about the role 
of (in this case SL(3)-) invariant polynomials associated with the system, the related question 
regarding Lagrangian structures (which exists for the all the lattice Gel'fand-Dikii hierarchy, 
cf. [E]), and the question about a possible existence of elliptic variants of the system. Since 
the lattice BSQ equations such as (II. ID . can sometimes exhibit interesting geometrical aspects, 
e.g. the emergence of the pentagram map, cf. [26] . matters regarding the discrete differential 
geometry in the sense of [9] may be worth considering as well. These and other issues will be 
relegated to future work. 
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